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Abstract 

We construct a one-parameter family of exact time-dependent solu- 
tions to 2+1 gravity with a negative cosmological constant and a mass- 
less minimally coupled scalar field as source. These solutions present 
a continuously self-similar (CSS) behavior near the central singularity, 
as observed in critical collapse, and an asymptotically AdS behavior 
at spatial infinity. We consider the linear perturbation analysis in this 
background, and discuss the crucial question of boundary conditions. 
These are tested in the special case where the scalar field decouples and 
the linear perturbations describe exactly the small-mass static BTZ 
black hole. In the case of genuine scalar perturbations, we find a grow- 
ing mode with a behavior characteristic of supercritical collapse, the 
spacelike singularity and apparent horizon appearing simultaneously 
and evolving towards the AdS boundary. Our boundary conditions 
lead to the value of the critical exponent 7 = 0.4. 
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1 Introduction 



The exploration of possible connections between seemingly different fields 
of physics has always attracted much interest and curiosity. In this respect, 
the observation first made by Choptuik jjj that gravitational collapse does 
exhibit many of the features typical of critical phenomena has opened a 
new of line of research (see for instance the review paper |§ and references 
therein). Specifically, Choptuik showed numerically, on the example of a 
gravitating spherically symmetric massless scalar field, that certain field 
configurations lying at the black hole threshold in the phase space of initial 
data for the gravity-matter system are characterized by universality, power- 
law scaling of the black hole mass and self-similarity (discrete self-similarity 
in this specific case). Recently, Pretorius and Choptuik Q have considered 
the same problem in a simplified context, namely the collapse of a massless 
and minimally coupled scalar field in 2+1 dimensional AdS spacetime. The 
inclusion of a negative cosmological constant is necessary for the existence of 
vacuum black hole solutions in 2+1 dimensions, the BTZ black holes By 
considering numerical collapsing configurations close to black hole formation, 
they obtained threshold solutions exhibiting characteristic critical features, 
namely power-law scaling and continuous self-similarity (CSS). Numerical 
calculations along similar lines were also performed by Husain and Olivier 

While the numerical investigation of critical collapse is by now a well- 
established branch of general relativity, far less is known on the purely an- 
alytical level. In 3+1 dimensions, a family of exact CSS solutions to the 
Einstein-massless scalar field model, first given by Roberts Q, was studied 
by various authors @-[10|. In particular, following the general analysis of 
linear perturbations of critical solutions JO, O], Frolov studied the linear 
perturbations of the threshold Roberts solution, and determined the mass- 
scaling exponent. In 2+1 dimensions, Garfinkle [13] found a class of CSS 
solutions to the gravitating massless scalar field model with vanishing cosmo- 
logical constant, and argued that one of these solutions should approximate 
the critical solution observed in || near the singularity. The drawback is 
that it is not possible to show that linear perturbations of such a solution 
lead to black hole formation, both because the solution of the perturba- 
tion problem necessitates boundary conditions which cannot be consistently 
stated for a solution which is known only in the vicinity of the singularity, 
and more fundamentally because black hole formation in 2+1 dimensions 
requires a cosmological constant, which does not appear in Garfinkle's CSS 
solutions. Subsequently, the present authors [Q found another CSS so- 
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lution to the same problem, and showed how to extend this solution to a 
quasi-CSS solution of the full problem with cosmological constant. It was 
then possible to perform the linear perturbation analysis of this extended 
solution, and to show that it did lead to black hole formation. 

In this paper we generalize the work of [ 14 1 to a new one-parameter class 



of CSS solutions, and significantly refine the analysis of black-hole forma- 
tion. In section 2 we review Garfinkle's one-parameter family of exact CSS 
solutions to the A = equations of motion. The question of extending these 
to solutions of the full A / system is addressed in section 3, where for 
a special parameter value we find an exact cosmological solution. Then, 
in sections 4 and 5, we use a limiting procedure to derive from Garfinkle's 
solutions a new class of A = CSS solutions, discuss their properties, and 
carry out their extension to solutions for the case with negative cosmological 
constant. The main contribution of this paper is represented by sections 6, 7 
and 8, where we consider the linear perturbation analysis in the background 
of these quasi-CSS solutions. The discussion of the most appropriate bound- 
ary conditions to be imposed on the perturbations, a separate analysis for 
the limiting case where the scalar field decouples, and a discussion of two 
possible scenarios for black hole formation in the linear approximation are 
then carried out. Finally, in section 9 we compare our results with those of 
previous numerical and analytical investigations. 

2 First class of CSS solutions 

The Einstein equations with a cosmological constant and a minimally cou- 
pled massless scalar field source are 

Gp, - Ag^ = kT^ , (2.1) 
where T^ v is the stress-energy tensor for the scalar field 

= d^dvt - ^d x (f>d x( f> , (2.2) 

the cosmological constant A = — l~ 2 is negative for asymptoticaly AdS space- 
time, and the gravitational constant k is assumed to be positive. 

The (2+l)-dimensional rotationally symmetric line element may be writ- 
ten in double null coordinates as: 

ds 2 = e 2a dudv -r 2 d8 2 , (2.3) 
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with metric functions a(u,v) and r(u,v). The Einstein equations ( |2.1| ) 
(where we take k = 1) and the associated scalar field equation are 
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For a given solution of these equations, the Ricci scalar is 

R = _6A + 4e- 2ff ^ J - 



(2.4) 

(2.5) 
(2.6) 
(2.7) 
(2.8) 



(2.9) 



Assuming A = 0, Garfinkle has found [13| the following family of exact 
CSS solutions to these equations 
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(2.10) 



depending on an arbitrary constant c and a scale ^4 > 0. These solutions are 
continuously self-similar with homothetic vector (u<9„ + vd{,). An equivalent 
form of these CSS solutions, obtained by making the transformation 



u 



(- U ) 2 1 , v = (y) 2q 



(l/2g = 1 - c 2 



(2.11) 



to the barred null coordinates (u,v), is 



ds 2 = -A{v q + {-u) q )^ 2q - l)/q dudv 
(j) = -2cln{v q + {-u) q ) . 



.(v 2 «-{-uf q fde\ 



(2.12) 



Garfinkle suggested that the line element ( [2.10 ) describes critical collapse 
in the sector r = — (u + v)/2 > 0, near the future point singularity u = v = 
0, where the effect of the cosmological constant can be neglected. The 
global spacetime structure of depends on the value of the parameter c 2 and 
is discussed in [14|. As shown in Jl3| ], the metric ( p. 12 ) can be extended 
through the surface v = only for q = n, where n is a positive integer. 

Let us mention here that these solutions belong to a larger class of solu- 
tions which may be generated from the Garfinkle solutions ( 2.1 0| ) or ( p. 12 ) by 



4 



the Ida-Morisawa transformation. Ida and Morisawa considered [15| three- 
dimensional Einstein-scalar fields with a hypersurface orthogonal spacelike 
Killing vector, parametrized by (|2.3D, and showed that the transformation 



(r', (/,(/)') = {r,a + b(j) + (6 2 /2) lnr, + Mnr) (2.13) 

generates a family (parametrized by b) of Einstein-scalar fields. Applied to 
( |2.12| ), this transformation leads to the solutions 

ds' 2 = -I'(^ + (_u)9)( 2c - fe ) 2 (-^ + (-u) 9 ) b2 dn^-i(^-(-n) 29 ) 2 ^ 2 , 
4' = (b-2c)ln(v q + (-u) q ) + bln(v cl -( y -u)< 1 ), (2.14) 

with q = 1/2(1 — c 2 ). In the case b = 2c, this solution (which was previously 
given in [JIJ], Eq. (2.16)) is again CSS. For 6/0, these solutions have a 
timelike central (r = 0) curvature singularity at v q — (—u) q = 0. They can 
again be extended through the surface v = for q = n. In the case q = n 
odd, the extension of a solution ( 2.14 ) with parameters (A, b, c) leads to a 
solution of the same class with parameters (A' = —A, b' = b — 2c, c' = — c). 



3 A class of quasi-CSS cosmologies 

For a consistent analytical treatment of critical collapse, we should be able 
to extend the A = CSS solutions ( |2~1~0| ) to A < quasi-CSS solutions, i.e. 
to construct (explicitly or implicitly) solutions of the full field equations with 
the behavior ( |2.10| ) near the singularity. However, we have been unable to 
find an ansatz which separates the variables and reduces the field equations 
to ordinary differential equations, except in the special case c 2 = 1/2 (q = 1). 
In this case the solution ( 2.12| ) with A = 1 may be written, after transforming 



to coordinates (R, T) defined by 

R-T R+T 



(3.1) 



V2 ' y/2 ' 

in the form of a FRW cosmology with flat spatial sections, 

ds 2 = T 2 (dT 2 - dR 2 - R 2 d9 2 ), 

$ = -v^lnT. (3.2) 

This form suggests searching for an extension by separating the equations 
in the variables R, T. With the metric parametrization 

ds 2 = e 2 *(dT 2 -dR 2 )-r 2 d8 2 , (3.3) 
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these equations read 
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with ' = d/&T and ' = d/dn. Let us make the ansatz 

(j> = A(T) + B(R), r = F(T)G(R). (3.5) 
The equations separate and lead to a FRW cosmology if 

G" = kG (3.6) 



(k constant). One finds that for the solution to reduce to ( |3,2| ) for A = one 
must have B constant, e.g. B = 0, the remaining equations being integrated 
by 



F 2 - a 2 12 - kF 2 - A(3 2 F 4 = 0, (3.7) 

where a and (3 are integration constants; without loss of generality we may 
choose a = \[2 (by rescaling times) and (5 = 1. The resulting metric 

ds 2 = F 2 (T){dT 2 - dR 2 - G 2 (R)d9 2 ), (3.8) 

where G(R) solves the equation G" = kG with the initial conditions G(0) = 
0, G"(0) = 1 for regularity, is of the FRW type, and reduces to ( |3.2| ) in the 
limit A — > if k — > in this limit. 

The discriminant (k 2 — 4A) of the effective potential in the last equation 
( |3.7|) is positive definite for A negative, so that this potential has two roots 
F 2 = T- < and F 2 = t+ > 0. It follows that this cosmology is time sym- 
metric around the turning point F 2 = r+, with initial and final singularities 
at F 2 = 0. Defining the new time variable r = F 2 (T), the line element ( |3.8| ) 
may be rewritten asQ 

ds2 = T7 TT V - < dR2 + G 2 (R)d9 2 ), (3.9) 

4 (t-t_) (t + - T) 

For G(i?) = const, (fc = 0, so that r_ = — r+), this solution reduces to the cosmological 
version of the static solution (4.4) with a — (put p_ = — p+ = — r+ and p 1 = t\ — r 2 ). 
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(A = -r 2 ) with 

G(R) = fi^sinhfiR (k = fi 2 >0), 

G(R) = R (k = 0), (3.10) 

G(R) = v^smvR (k = -v 2 < 0). 

In the special case k = 0, this line element may be rewritten as 

ds 2 = F 2 L^jp ~ dR 2 - R 2 d6 2 ^ , (3.11) 
showing clearly the restoration of self-similarity in the limit I — * oo. 

4 Second class of CSS solutions 



From the Garfinkle class of CSS solutions (2.10), we may derive a new class 
of CSS solutions by a generalization of the construction of [14]. The double 
null ansatz fl2.3|) is invariant under boosts u — > au, v — > a~ x v. Rescaling 
angles by 9 — > (2/a)6 and putting A = a~ 2 ° 2 A, we obtain from (2.10), in 
the limit of an infinite boost a — > oo, 

ds 2 = -A{^fdudv-u 2 de 2 , 

v 

ct> = -cln(-n). (4.1) 

These new solutions are obviously, as the original Garfinkle solutions, con- 
tinuously self-similar with homothetic vector (u<9„ + vd{,). A simpler form 
of these CSS solutions may be obtained by transforming to new null coor- 
dinates 

u = -(-u) l+c \ v = v 1 ~ c \ (4.2) 
which yields (with the normalization A = c 4 — 1) 

ds 2 = dudv-(-u) 2 ^ 1+c ^d9 2 , 

<t> = --^\n{-u). (4.3) 
I + c z 

For c 2 = 1 we recover the "new CSS solution" of [14]. 



Interestingly, the CSS solutions ( |4.l| ) may also be obtained by submit- 
ting the singular static solutions of the cosmological Einstein-scalar field 
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equations to an infinite boost near the singularity. These solutions, derived 
in H, 

rfs 2 = A\p-p + \V 2 + a \p-p„\V 2 - a dt 2 - ^-\p- P+ \V 2 - a \p-p„\V 2+a d9 2 



dp 2 ll -4a 2 



+ 4A0 ,-;; )(p - p - r *-v^nraJ (4 - 4) 

(with —1/2 < a < 1/2, and p± = 0(A -1 )) are singular for p = p±. In the 
neighborhood of one of these singularities, e.g. /)+, the transformation p = 
p + + a; 4 /*- 1-2 ") leads (after appropriate rescalings) for x — > 0, or equivalently 
for A — > 0, to the exact A = solution [fLTl] 



dfi 2 = x 2(l+2a)/(l-2a) ({ft 2 _ ^2) _ ^2^ = ^ 



2(1 + 2a) 

-lnx. (4.5) 

2a 



Transforming to null coordinates u = t — x, v = t + x, and boosting this 
solution by u — » au, {) — ► a" 1 ^, we obtain in the limit of an infinite boost 
(again after appropriate rescalings) 



ds > = (-uf(^/^) d udv-u 2 de 2 , cf> = ^ffi±Ml n |u|. (4.6) 



After a suitable transformation on the coordinate v, this is found to coincide 
with the CSS solution ( |4,1| ) provided we identify 

9 1 + 2a . . 

For a = —1/2 (+1/2) the solution ( [4.4| ) is no longer singular for p = p + 
(p_) but reduces to the BTZ black hole 

ds 2 = {r 2 ,l 2 - M)dt 2 - {r2/ f_ - - r 2 d9 2 , (4.8) 

with cj) = 0. Near r = 0, this reduces for M > to the Minkowski cosmology 
ds 2 = M _1 dr 2 — Mdt 2 — r 2 d6 2 , which after an infinite boost leads to the 
metric rt4.3|) with c = 0, 



ds 2 = dudv-u 2 d0 2 . (4.9) 

The infinite boost limit corresponds physically to the massless limit, and so 
one would expect that the metric (O) can alternatively be derived from 
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(4.8) by taking the limit M — > (this time without an infinite boost). 
Indeed, the extreme (M = 0) BTZ black hole, or "vacuum solution" 



ds 2 = (r 2 /l 2 )dt 2 - 4^777 - r 2 d6 2 , ( 1. 10) 



may be rewritten (with x = l 2 /r) as 

I 2 



ds 2 = ^(dt 2 -dx 2 -l 2 d6 2 
x z 

Al 2 



{U _ r)2 -(dudv-i 2 de 2 ) (4.ii) 

in null coordinates U = t — x, V = t + x. The transformation 

U = -2l 2 /u, V = v/2 (4.12) 



transforms ( 4.11 ) to 

ds 2 = ——^(dudv - u 2 d6 2 ). (4.13) 

(1 + uv/Al 2 ) 2 

This goes over to the c 2 = CSS solution (|4.9|) in the limit A — > (/ — > cxd) 
or, equivalently, in the near-horizon limit r — > (n — > 0). 

The c 2 = spacetime (4.9) is Ricci-flat, and so (in 2+1 dimensions) is 
flat. Therefore it is locally equivalent to Minkowski spacetime. The explicit 
local transformation from the Minkowski metric written in the double null 
form 

ds 2 = dudw -dy 2 (4.14) 



to the metric (4.9) is 

u = u, y = u9, w = v + u9 2 . (4-15) 

The rotationally symmetric c 2 = CSS spacetime is obtained by periodically 
identifying 8 with 0+2tt. From this vacuum solution (r = —u, a = 0, <j> = 0), 
the entire second class (4.1) of c 2 / CSS solutions may be generated by 
the Ida-Morisawa transformation fl2.13j ) with b = c. 

For all c / 0, the CSS solutions ([D^) have a high degree of symmetry, 
with the four Killing vectors 

1 

L\ = vd v - ud u + 9 0d e , 
& + 1 

£ 2 = -2^^ed v + (-u)^d e , 

& + 1 

L 4 = 3 e . (4.16) 
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For c 2 = 1, the above form of L 2 is replaced by 



L 2 = 9d v -Un{-u)d e . (4.17) 

As in the case c 2 = 1 fli"4|| , the Lie algebra generated by the Killing vectors 
( |4.16| ) is solvable. In the case c 2 = 0, the flat metric (|4.9| ) admits of course 
six Killing vectors generating the (2+l)-dimensional Poincare algebra. 

Again as in the case c 2 = 1, the metric (|4.3|) is devoid of scalar curvature 
singularities. The nature of the corresponding geometry may be understood 
from the study of geodesic motion in this spacetime. The geodesic equations 
are integrated by 

u = tt, (-u) 2/ ( c2+1) e = -I, TTV + W = e, (4.18) 

where ir, I and e are constants of the motion. The spacetime is extendible 
to a geodesically complete spacetime in two cases: 

a) c 2 = 0. This is the near-horizon limit (|4.9| ) of the BTZ vacuum ( 4,13| ), 
and is diffeomorphic to the Minkowski metric, as explained above. 

b) c 2 — > 00. In this limit the scalar field decouples and the metric 
reduces to the geodesically complete cylindrical Minkowski metric. 

In the other cases, the third equation ( 4.18| ) integrates to 



v = ~ l Y 2 +1 h -u) (c2 - mc2+1) + ™nst. , (4.19) 

TT Z 7T Z (C Z — 1) 

showing that for c 2 < 1 nonradial geodesies terminate at u = 0,v — > +00, 
while radial geodesies (I = 0) can be continued through the null line u = 0, 
which they cross at finite v, to u — > +00, only the endpoint v — > +00 of the 
null line u = being singular. In the case c 2 = 1 (treated in 0]), ( [4.18| ) is 
replaced by 

£ I 2 

v = — -u r ln(— u) + const. , (4.20) 

leading to the same conclusion as for c 2 < 1. On the other hand, for c 2 > 1 
nonradial geodesies terminate at u = 0, v = const., and the whole line 
u = corresponds to a singularity of the geometry. So the c 2 = 1 solution 
is an extreme solution in the manifold of CSS solutions (|4.3| ), lying at the 
threshold between two classes of solutions (0 < c 2 < 1 and c 2 > 1) differing 
by their global spacetime geometry (Fig. 1). 
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5 Quasi-CSS extensions of the second class 

In this section we proceed to extend the second class of A = CSS solutions 
fl4.3|) to exact solutions of the full A < equations (the self-similarity being 



then broken by the cosmological constant). Following [14] we make the 
ansatz 

ds 2 = e 2u{x) dudv- (-u)^p 2 (x)d9 2 , (f> = — JL-hx\ u \+ib(x) , (5.1) 

c z + 1 

with x = uv. This reduces to the CSS solution ( p~3| ) for p = 1, v = ip = 0, 
and preserves the Killing subalgebra (L\,L4). Inserting this ansatz into the 
field equations (|2.4|)-(|2~8|) leads to the system 

V' + ^P' = (5.2) 

2(xz/' + z/)+^W- Vr?) = T^' ( 5 - 3 ) 

c z + 1 2 

x 2 (-p" + 2pV - p ^) + -?— x (-pi + p (y> + «//)) = (5.4) 



c 2 + l 



p" + 2p'v' - pip' 2 = (5.5) 



2x{piP>)> + 2 4±lp^' = ^— //(5.6) 
c z + 1 <r + 1 



(' = d/dx). A simple first integral, obtained by comparing ( |5.4| ) and ( |5.5| ), 
is 

p = e^^ (5.7) 



(with the integration constant set to 1 by the boundary conditions (5.8)). 
The boundary conditions 

p(0) = l, K0)=0, V(0) = (5.8) 

lead to a unique solution reducing to Q4.3; ) near u = 0. The small x behavior 
of this solution is 

c 2 + 1 



1 H r-s rAx, 

2(c 2 + 2) 

(c 2 + l)(c 2 + 3) 

2(c 2 + 2)(2c 2 + 3) 



Ax, (5.9) 



^ ~ 2( C 2 + 2)(2 C 2 + 3) AX • 
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The analysis of [14] can be straightforwardly extended here, leading to 
the conclusion that when x decreases from x = 0, the functions p and e 2u 
increase indefinitely, going to infinity for a finite value x = x\. The numerical 



solution of the system (|5.2|) -( |oT6D with boundary conditions fl5.8|) leads to 
values xi(c 2 ) equal to — 4Z 2 for c 2 = and for c 2 — > oo, and remaining of 
the order of — 4/ 2 for finite c 2 , the maximum value being xi(l) = — 3.876Z 2 . 
The behavior of the metric functions near x\ is found to be 

(\ c 2 c 2 xln(x) \ 

9 = Pl \~x~ + 2(c 2 + l)xi ~ 12(c 2 + l) 2 x 2 + "') 

2 „ _ 4xi ( c 2 x 2 ln(x) \ 

6 ~ Ax 2 l v 1+ 12(c 2 + l) 2 x 2+ - y l (5 ' 1Uj 

cx cx 2 , , . 

^ = ^1 + n( 2 i 1\ / 2 , 1N2 2 m ( X ) + ■•• 

2(c^ + ljxi 8(c^ + lj^xf 
(x = x — xi). Changing to coordinates (T,R) defined by 

u = _^+l e R-T^ v = l l-^ e R+T^ (5n) 

we obtain from ( |5.10| ) the leading asymptotic behaviors near R — R\ 
ds 2 ~ ^ -, 2 (df 2 - cf£ 2 - e^ (?1 - f W) , = 0i + cT/(c 2 + 1) 

(Ri - Ry 

(5.12) 

(R — R\ ~ x/2xi). Making the further coordinate transformation, 

fl-fli = -(c 2 + l)/XT, T-Ti = (c 2 + l)ln(T/(c 2 + l)), (5.13) 
we arrive at the asymptotic form 

ds 2 ~l 2 (^X 2 dT 2 - < ^-X 2 de 2 ^j, = 1 + clnT, (5.14) 

showing that the metric is asymptotically AdS (with logarithmic subdomi- 
nant terms). The corresponding conformal diagrams for the cases < c 2 < 1 
and c 2 > 1 are shown in Fig. 2. 

For c 2 = or c 2 ^ oo, the solution of Eq. ( pjq ) is pi\)' = ax p , with a 
constant and p = 3/2 for c 2 = or p = 1 for c 2 — > oo, so that the boundary 
conditions ( |5.8| ) enforce a = 0, i.e. '0 = 0, and the extended metric (5.1) is 
a vacuum metric. For c 2 = this metric is 

ds 2 = -. , , (dudv — u 2 d6 2 ), (5.15) 

(l + uv/Al 2 ) 2 



12 



which we recognize as the BTZ vacuum ( f4.13| ). For c 2 — ► oo, the extended 



metric 2 

ds 2 = —^(dudv - (1 - uv/4l 2 ) 2 d6 2 ), (5.16) 

(1 + uv/Al 2 ) 2 

is obtained from the BTZ metric Q4.8| ) by transforming first to coordinates 
(R, T) with 

r = M 1/2 l coth(M 1/2 i?//), t = T, (5.17) 

which leads to 

^ = S inh 2 (M^R/l) (dT2 " ^ ~ M ^ coth2 ( Ml/2i? / / ) ( 5 - 18 ) 

then, choosing the convenient normalization Ml 2 = 1, to coordinates (u, v) 
defined by 

u = -2le^l l \ v = 2le l - R+T ^ l \ (5.19) 



So the family of time-dependent quasi-CSS solutions ( |5.1| ) interpolates be- 
tween the BTZ vacuum for c 2 = and the massive BTZ metric for c 2 — > oo. 



6 Perturbations: general setup 



Now we study linear perturbations of the quasi-CSS solutions (5.1) along the 



lines of the analysis of |9|, |l0|, 14 1 . The relevant parameter in critical collapse 



being the retarded time |l3| T = — ln(— u) = — (l/(c 2 + l))ln(— u), we 
expand these perturbations in modes proportional to e kT = {—u)~ k ^ c +l \ 
with k a complex constant. Only those modes with Re(k) > 0, which grow 
when T — * +oo, will possibly lead to black hole formation. Keeping only 
one mode, we decompose the perturbed fields as 

r = {-u) 1 ^ c2+1 \p{x) + {-u)- k ^ c2+1 ^{x)), 

\n\u\+^{x) + (-u)- k/{c2+1) 4>(x), (6.1) 



c 2 + l 

a = v {x) + (-n)- fe /( c2+1 )<?(x) 



The linearization of the Einstein equations Q2.4| )-( |2T8| ) in the perturbations 
r, 4>, a leads to the system 



xr H 5 — r =— e (r + 2/ocr), 

cr + 1 



2 In this case the derivation of ( p.7| ) breaks down. The solution ( p.lq) is obtained by first 
^r 



i nteg rating; ( p.5| ) to p' = const, e then integrating ( p.2| ) with the boundary conditions 
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c 2 + 1 J Y C 2 + 1 
f + pxa' — k I xp' + 



2x a" + 2(C ' 9 + 1 fc) a' = ke 2v ~a - ( 2x^ - 
c z + 1 V 

-(1 - k)xr> + ((1 - k)xv> + 

— p (cx(f> — k ( 2 — - — x0'^ — cxip'r 
2(p'a + u'f') - r" = ip'(2 P 4>' + if>'r), 



+ 



c 2 + l 



, 2c 2 + 3 - 2/c 



c 2 + l 



c 2 + 1 



p'^ + ( 2xip' 



c 2 + l 



/ „ 2c 2 + 3 -k A 
+ ( 2x^" + c2 + 1 V J f = 0. 



This linear differential system is of order four (Eqs. (6.4) and ( |6.5[ ) are 
constraints, while the perturbed scalar field equation ( |6.6D is a consequence 
of the perturbed Einstein equations). An exact solution of the system (|6.2j)- 
(|6 .6|) is given by 



f fc (x) 


= a(— 


1+k/{c 2 + 1)p , {x) ^ 




= a(— x] 


i+*/( C 2 +iy (x ) , 


a k (x) 


= <*((- 


x )l+V(c 2 +l)^ (a . 



(6.7) 



c 2 + 1 + k 
2(c 2 + 1) ' 



xfc/(c 2 +l) 



The corresponding first-order perturbed fields (|6.1| ) are generated from the 
unperturbed fields by the gauge transformation v — > u-qu 1+<: ^ c +1 \ So, up 
to gauge transformations, the general solution of this system depends only 
on three integration constants, which should be determined by enforcing 
appropriate boundary conditions. 

At this point we should emphasize that the question of which boundary 
conditions are truly "appropriate" is rather subtle in the general-relativistic 
context (for instance, in the related problem of scalar field collapse in 3+1 
dimensions, different boundary conditions have been chosen in || and JlOjl), 
yet crucial for a proper determination of the eigenmodes. 

The three boundaries of our quasi-CSS solutions are u = 0, v = (or x = 
0), and x = xi (the AdS boundary). At the boundary u = 0, corresponding 
to the original centre of the unperturbed quasi-CSS solution, it seems natural 
to impose that (A) the perturbed metric component r does not diverge too 



quickly, which would conflict with the linear approximation [10|. However 
this is subject to some ambiguity. For instance, should one consider the 
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limit u — ► with v held fixed, or rather with x = uv held fixed? Also, the 
perturbed spacetime may develop an apparent horizon hiding the singularity 
at u = 0. In such a case it may seem preferable to view this apparent 
horizon, rather than the line u = 0, as the natural boundary of the perturbed 
spacetime on which to impose a boundary condition. 

On the second boundary v = 0, it seems natural to require that (Bl) the 
perturbed solution matches smoothly the original quasi-CSS solution ||] 

f(0) = 0, cr(0) = 0, 0(0) = (6.8) 

(on account of ( |5.8| ), the perturbed solution then also matches smoothly 
the original CSS solution). However conditions (pjjj) do not ensure that 
the curvature tensor remains finite on u = 0, which seems essential for a 
smooth matching. This last requirement is automatically satisfied if the 
"weak" conditions ( |6.8| ) are supplemented with the the "strong" condition 
(B2) that the perturbations be (as the unperturbed quasi-CSS solution) 



analytic in v [14]. 



Finally, we shall require that (C) the perturbations grow slowly (i.e. not 
faster than the original fields) as the third boundary x = x% is approached 
U , this condition ensuring the validity of the linear approximation near the 
AdS boundary. Note that, owing to the gauge freedom, it is only necessary 
that these various boundary conditions be satisfied up to gauge transforma- 
tions, i.e. up to the addition of a gauge perturbation (|6.7D . 

To enforce boundary conditions on the boundary v = 0, we must first 
investigate the behavior of the solutions of the linearized system near the 
boundary x = 0. Generalizing the analysis of ]l4j| , we assume the power-law 
behavior 

r(x) oc {-xf (6.9) 
(p constant). For the purpose of the determination of the exponent p, Eqs. 



fl6.2|), (|6.3D and (6.5) can be approximated near x = as 



_„ , c 2 + 2 - k 

c 2 + 1 
,„ t c 2 + 1 - k 

c 2 + l ° 2(<-' + i) 



xf" ^ -I f'~Aor, (6.10) 



g^A^-f^- ; 9 C(C ' + ^ rA0' (6.12) 

c 2 + 2 (c 2 + 2)(2c 2 + 3) V K ' 

(in ( |6.12| ) we have replaced the unperturbed fields by their small x behaviors 
( |5.9|) ). The discussion of this system is simpler in the special (vacuum) case 
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c 2 = 0, which we shall consider in the next section, leaving the discussion of 
the general (scalar) case c 2 7^ to Sect. 8. 



7 BTZ black holes as vacuum perturbations 

In the case c 2 = 0, we have seen that the quasi-CSS solution is the BTZ 



vacuum, with ip = 0. Then, the linearized equations (|6.2| )- (|6.6|) decouple 
into the linearized sourceless Einstein equations and the linearized scalar 
field equation on the BTZ vacuum background. The solution of the three- 
dimensional cosmological Einstein equations being essentially unique, up to 
diffeomorphisms, the three-parameter family of perturbations of the BTZ 
vacuum must include the small-mass BTZ black holes, with the mass as 
perturbation parameter (amplitude), as we now check explicitly. 

When c 2 = 0, Eqs. (|6.1l| ) and (|6.12| ) decouple and must be satisfied 
separately. Then a can be eliminated between ( 6.10| ) and ( |6.12| ), leading to 
the third order equation 

xr'" + (1 - k)f" ~ 0, (7.1) 
which implies the secular equation 

p(p- l)(p- 1 - k) = 0. (7.2) 



We note that the root p = 1+k corresponds to the gauge mode (|6.7| ). Dis- 
carding this mode, we obtain from the full equations ( |6.2j )-( |oT6|) the behavior 
of the perturbations near x = in terms of two integration constants A, B: 

r(x) ~A + B(-x), (7.3) 
~ -J- (2-k)A~ 1 B + 0(x). (7.4) 



This can only satisfy the boundary conditions (p.q) at x = if 

A = 0, k = 2. (7.5) 



It follows that in this case the solution of the linearized system (|6.2|) -( |6T5j ) 
with the initial conditions fl6.8| ) is, up to a gauge transformation, unique. 
We now show that this solution is simply a linearization of the BTZ metric. 



The BTZ metric ( 5.18; ) may be written in double-null form as 



4/ 2 dU dV 



ds z = 27 - Ml z coth z {U -V)dd\ (7.6) 

sinh (U — V) 
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with U = n(T + R), V = p(T - R), and fx = M 1 / 2 /2l. This may again be 
rewritten, in terms of the new null coordinates 



2^ 2 cothC7, v = (2///) tanh V, (7.7) 



as 



ds 2 = Sinh i 2 2 ^ C r ° Sh l y du dv - 4/i 2 ^ coth 2 (t/ - V) dB 2 



s\nh.\U - V) 



' (dudv - (u + (M/A)v) 2 d6 2 ). (7.1 



(l + uv/Al 2 ) 2 



This form of the BTZ metric again reduces to the vacuum CSS metric 
in the infinite boost limit (which amounts to taking v — > 0) and, furthermore, 
can obviously be considered for small masses as a perturbation of the BTZ 
vacuum ( 4.13| ). From ( [7,8| ) the metric function r is 



r = p (-u - (M/A)v) 

= -u{p+{M/4)u- 2 xp), (7.9) 

with x = uv, and p(x) = (1 + x/4Z 2 ) -1 . This is of the form of ( |6.1| ) for the 
eigenmode k = 2, as found above. On the other hand, the metric function 
e 2u in (|7.8| ) does not depend on the perturbation parameter M/4, so that 
the perturbations f and a are simply given by 

M 

r = —xp, a = 0, (7.10) 

which satisfy the condition of slow growth at the AdS boundary x — ► — 4l 2 
(where f ~ —Ml 2 p), as well as the requirement of finiteness of r for u — » 
with u held fixed (but not for u — > with x held fixed), and are easily 
checked to solve exactly the vacuum system (|6.2D - (|6.5[ ). 

It is instructive to analyze more closely in this (obviously very special) 
case the phenomenon of black hole formation. First we note that the un- 
perturbed metric ( 4.13f) can be extended beyond u = — oo by the coordinate 



transformation u = —u , leading to the metric 

ds<2 = 7 l -TJTTv? ( d udv - d9 2 ), (7.11) 

[u — v/4l z ) z 

which is regular in the triangle u — v/Al 2 > 0, the other boundaries of 
this triangle being the apparent horizons r „ = (u = 0) and r u = 
(v = — oo) through which the metric can be periodically extended to the 
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geodesically complete BTZ vacuum spacetime (Fig. 3). Bearing in mind 



that the perturbed metric (7J)) can similarly be extended beyond u = — oo 
and v = — oo, and is obviously extendible to u > 0, we find that the effect 
of the perturbation is twofold (Fig. 4). First, the two apparent null central 
singularities u = and v = — oo are replaced by two genuine spacelike 
central causal singularities @ u+(M/A)v = 0. Second, the two null apparent 
horizons, displaced to 

r, v = -u\p> + (M/4) U - 2 (x P y)=p\^- T )=0, 

Mv 2 

r, u = -(p + xp' + {M/4:)u- 2 (-xp + x(xpy)) = -p 2 (l-— F )=0, 



161 2 



(7.12) 



effectively shield the final and initial singularities. 

While the situation here is static, rather than dynamical, this formation 
of apparent horizons in the linearized approximation (here exact) is formally 
similar to that observed in critical collapse. Recall that near-critical collapse 
is characterized by a critical exponent 7, defined by the scaling relation 

Qoc\p-p*r, (7.13) 

for a quantity Q with dimension s depending on a parameter p (with p = p* 
for the critical solution). In the present case, it has been shown |l8[| that the 
BTZ vacuum is a critical point of the spinless BTZ black hole, and that the 
phase transition from extremal (M = 0) to nonextremal BTZ black holes is 
second order, so that a similar scaling relation applies. Choosing for Q the 
radius 

r AH = M x l 2 l (7.14) 

of the apparent horizon (s = 1), and identifying p — p* with the perturbation 
amplitude M/4, we obtain from ( [7. 14 ) the value of the BTZ critical exponent 



7 = 1/2, (7.15) 

in agreement with the value derived previously either from the dynamical 
analysis of black- hole formation in 2+1 dimensions from spherical dust col- 
lapse [19[ or point particle collisions [20], or from the consideration of phase 



transitions from non-extremal to extremal black holes IllHl or from static to 



rotating black holes [21]. 
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8 Scalar perturbations and black hole formation 



Let us return to the case c 2 7^ where the scalar field does not decouple. 
Eliminating now the functions a and between the three equations ( 6.10| )- 
( |6.12| ), we arrive at the fourth-order equation 

9 .... 2k -3c 2 -7/2 ... (k - c 2 - l)(k - c 2 - 3/2) „ , , 

x f 2 1 1 x ? + 7 2 , ^2 ~ f - °> S- 1 

c z + 1 (c z + l) z 

which leads to the secular equation 

x f k \ f k + c 2 + l/2\ , N 



Again, the root p = 1 + k/(c 2 + 1) corresponds to the gauge mode (6.7). 
Discarding this mode, we obtain the general behavior of the perturbations 
near x = in terms of three integration constants A, B, C: 



d{x) \ \ A-'B + 0(x) 



f(x) ~ A + B(-x) + AC(-x) , (8.3) 

,4 c 2 + 2 - fc 

- (c2+;/ ^ + 1/2 ' c '-»^- 

~. , {2-k-c 2 ) A c 2 + 2-fc , 

4>(x) ~ ^ ^ + c(c2 - A x £? + O(x) 

,(c 2 + 3/2)(. + C 2 + l/2) c( _ x) ^ ) 



c(c 2 + l) 2 

fc + e 2 + l/2 

for fe 7^ 1/2 (in the degenerate case k = 1/2, the power law (— x) c ' 2 + 1 is 
replaced by a logarithm). 

From these behaviors, we find that the smooth matching conditions (Bl) 
at v = are satisfied if either 

A = 0, B = (Re(k) > 1/2), (8.6) 



or 



,4 = 0, k = c 2 + 2. (8.7) 



The stronger condition (B2) of analyticity in v then implies, in the first 
case, that k = n(c 2 + 1) + 1/2, where n is a positive integer, and, in the 
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second case, that C = 0. So the boundary conditions (Bl) and (B2) have 
two solutions: 



a) A = B = 0, k = n(c 2 + 1) + 1/2, (8.8) 

b) A = C = 0, k = c 2 + 2. (8.9) 

When u —* with f held fixed, the corresponding metric perturbations 
(— -u)( 1-fc )/( c +1 )f go to zero (in case a)) or to a finite limit (in case b)), so 
that in all cases the perturbed r has a finite limit and the boundary condition 
(A) is automatically satisfied. If on the other hand one takes the limit u — > 
with x held fixed, then in all cases r diverges as (— ( c +1 \ In case b) 
the divergence is the same (of the order (— -u) -1 ) as in the BTZ case, while 
it is weaker in case a) with n = 1, and stronger in case a) with n > 2. So 
a boundary condition that r does not diverge more than in the BTZ case 
when u — ► with x held fixed would select the eigenvalues k = c 2 + 3/2 
(case a) with n = 1), or k = c 2 + 2 ( case b)). 

Two of the three integration constants A, B, C being now fixed, the first- 
order perturbation is now completely fixed up to an arbitrary scale, and (as 



discussed after (6.7)) within a gauge transformation. So we now must check 
that the natural condition (C) at the AdS boundary is indeed satisfied up 
to a gauge transformation. The leading behaviour of the background at the 



AdS boundary (a; — > x\) is, from Eqs. ( |5.10 ), 



P^-^, e 2i W^) ? 1 ^Vi- (8.10) 
x — x\ \ A / (x — x\Y 

We again assume a power-law behavior 

d~bx q (8.11) 



[x = x — x\). Then Eq. ( [6.3] ), where </> can in first approximation be 
neglected, gives 

<z(g-l) = 2, (8.12) 
i.e. q = — 1 or q = 2. Accordingly, Eq. ( |6.2j ) reduces near x = to 

f" - 2x~ 2 f ~ 46pix 9 " 3 . (8.13) 

If q = —1, the behavior of the solution is governed by the right-hand side, 



i.e. 

-2 



r ~ (8.14) 
which apparently violates the boundary condition (C) for x — > x±. However 



the behaviors ( p. 11 ) with q = — 1 and ( p. 14 ) are precisely those of the gauge 
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mode ( |6.7[ ) near the AdS boundary, so that they can be transformed away 
by a gauge transformation, i.e. we can find a gauged] in which the behavior 
( |8.14| ) is replaced by 

r~-^-, (8.15) 
x — X\ 

(E constant) which is consistent with the boundary condition (C), and is an 
admissible small perturbation if its amplitude is small enough, J? C 1 A 
finer analysis, where cf> is not neglected in (|6.3| ), then leads to the behavior 
a oc x 2 In \x\, and to the behavior ( |3.15|) for f. 

Black hole formation is characterized by the appearance of a central 
singularity hidden by an apparent horizon. The definition of these notions 
can be quite ambiguous in linearized perturbation theory. For our present 
purpose we will consider, rather than the centre r = 0, the coordinate sin- 
gularity of the perturbed metric 

1 k 

^=e 2a r = (-u)^e Mx) [p{x) + (-uT^T (r(x) + 2p(x)a(x))] = 0. 

(8.16) 

Let us discuss the solution of this equation near x = 0, considering separately 
the two cases a) and b). In case a) (A = B = 0), we see from Eqs. ( |8.3| ) and 
( |8.4|) that f + 2/9CT is dominated by 2a, so that the coordinate singularity 
develops for 

^--V ^ + W/^ + W cw'SS, (8.17) 
which may be rewritten as 

i 



(-u)*&+v ocCv n . (8.18) 

For the lowest value n = 1 (k = c 2 + 3/2), we see that if C > there is (to 
the order considered) no coordinate singularity for v < 0, while a spacelike 
coordinate singularity develops for v > 0. Note that the consideration of 
the centre r = would have given quite different results, namely an eternal 

k 

timelike centre at (—u) c ^+ 1 = — r = —ACx (for n = 1). The location of 
the trapping horizon is given by 

<P" -j-2 k 

r }V = (-u)^+i(p' + (- u y^r') = 0, (8.19) 



3 This gauge transformation is equivalent to a shift in the position x\ of the AdS bound- 
ary. 
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or 

( _ )3fe ^ 2 (^ + 2K fc++ ^ + l/2) c( _ x) ^ (8 20) 

near x = 0. For n = 1, C > 0, this spacelike trapping horizon develops for 
v > 0, just as the spacelike coordinate singularity ( |3,17 ) which it shields for 



all c 2 (Fig. 5). While we are unable to show the formation of a curvature 
singularity without a nonperturbative analysis, it is suggestive that this 
formation of a shielded coordinate singularity is precisely what would be 
expected for a genuine curvature singularity. As the AdS boundary x = x\ 
is approached, it follows from the asymptotic behavior f ~ Ep + F that 
both the coordinate singularity and the trapping horizon approach the AdS 
boundary for a constant value of u, with a finite limit for the radius of the 
apparent horizon 



p \ ( ,r\ ( F 



-)E l l k . (8.21) 



Note that it is essential for this trapping horizon to form that the cosmo- 



logical constant be non-zero. For A = 0, Eq. ( 8.17 ) for the coordinate 
singularity becomes exact, while on the other hand p = 1 and f = every- 
where so that r jV is identically zero. 

In case b) (A = C = 0, k = c 2 +2), a finer analysis leads to the behaviours 
(consistent with the exact solution ( |7 . 1 0| ) ) near x = 

f ~ -B(x + j(l + c 2 )x 2 + ■■■), (8.22) 

c 2 

or ~ ~B(-x + •••), (8.23) 

~ £(<» + •••). (8.24) 

From these, it follows that that if B < 0, there is no coordinate singularity 
for v < 0, a spacelike coordinate singularity developing for v > according 
to 

(_ U )7^TI ~ + c 2) w . ( 8 .25) 

On the other hand, the location of the apparent horizon is given near x = 

(-u)P^ =~~B(l + ( 7 -^)s), (8.26) 
P 

with S = 2£(c 2 + 2)/A(c 2 + l) > 0, /? = -A(c 2 + l)/2, 7 = -A(c 2 + l)(4c 2 + 
9)/2(2c 2 + 3) 2 , < 7 < (3 (the first order behaviour of p' has been computed 
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from Eq. ( |5.2| )). This apparent horizon preexists the singularity (Fig. 6), 
its radius varying as 

r AH ~ [5(1 + (7 - 0)xp k (l + 0(x 2 )), (8.27) 

i.e. increasing with decreasing x, as suits a collapsing configuration. The 
behaviour of the apparent horizon in the region x > is less clear. In 
particular, it would seem to be asymptotic to the null line u = for some 
finite value x = xq (corresponding to v — ► — oo) , but at this point tah would 
diverge, in contradiction with our condition (A) for the validity of the linear 
perturbation analysis. A more refined calculation, including higher order 
perturbations, would be needed in order to clarify this issue. 

Now we inquire whether the apparent horizon of the perturbed spacetime 
is regular. The non-constant part of the Ricci scalar ( |2.9| ) is, to first order, 
the product of three factors, 

R + 6A = Ae- 2a (/) tU ()) :V = 4e~ 2u [1 - 2(-u)~^W] • 

+ xi// + (-tt)~^i( 9 fc i + xfi)}- 



c 2 + l c 2 + l 

[if/+ {-uf^#\. (8.28) 

Possible curvature singularities for v = having being excluded by the an- 
alyticity condition (B2), perturbative curvature singularities can only arise 
from the divergence of any factor of ( |3.28| ) on the original null singularity 
u = 0. Leaving aside case b), where as mentioned above it is not clear 
whether the apparent horizon intersects the line u = 0, we concentrate on 
case a). The perturbative Ricci scalar ( 8.28] ) evaluated on the apparent 
horizon ( |8.19| ) is 



+ xtf - // ?+1 V ~ ] W ~ ■ (8.29) 



+ xd>' 



c 2 + l 



The last factor <j> v diverges as x -1 when x — > 0, while the other factors go 
to finite limits, so that the apparent horizon will generically be singular at 
its onset x = 0, unless the product e~ 2a 4> u accidentally vanishes, 



e~ 2a (/>,u\AH(x = 0) = 0. (8.30) 



4 This again is due to a deficiency of the perturbative approach: the perturbative co- 
ordinate singularity due essentially to the vanishing of e 2<T may, as mentioned above, be a 
more correct indicator of a curvature singularity. 
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We suggest that this condition of "regularity at birth" of the appar- 
ent horizon is the crucial boundary condition (D) which will determine the 
eigenvalue of k. Specifically, for x — > 0, 

/ c 2 + 3/2\/ (2A;-l/2)(c 2 + 3/2)\ , 
flAg + 6Aoc (l + ^f j (c- 2 / ( J' +2) > J,-. (8.31) 

The first factor is positive, while the second factor vanishes for the value 

*° = i + c ?TW (8 ' 32) 

Unfortunately, it is easy to see that ko < c 2 + 3/2 for any c, so that the 
analyticity condition (B2) is not satisfied and consequently the Ricci scalar 
diverges on the line v = 0. Nevertheless we believe that the mechanism 
proposed here for the determination of k is basically correct, and that it 
is its implementation in the linear approximation which is responsible for 
the problem just mentioned. Indeed, a way to overcome this problem is 
to play with the ambiguities of the linear approximation and linearize the 
product e _2cr jU before evaluating it on the apparent horizon. Then, ( |3.31| ) 
is replaced by 

B _i_ K A ( i c ( c2 + 3 / 2 ) (2fc-l/2)(c 2 + 3/2) ^ _! 

Rah + 6A<x\c+— jt ^ 2^Pl) ) X ' (8 - 33) 

which vanishes for the value 

c 2 + 3/2 y ' 



This is smaller than (2c 2 +5/2) (the value ( |8.8| ) for n = 2), while the equation 

kx = c 2 + 3/2 (8.35) 

(n = 1) has only one real solution c 2 ~ 1.045, leading to k ~ 2.545. While 
the linearization trick we have used to obtain this value is hard to justify, it 
is nevertheless suggestive that we have found a set of boundary conditions 
(Bl), (B2), and (D) which determine a single value for the couple (c, k), i.e. 
a unique critical solution with a unique growing mode, and that this value 
c 2 ~ 1.045 is very close to the value c 2 = 1 suggested as a critical value by 
the global analysis of our family of CSS solutions (Figs. 1 and 2). Let us 
also recall that Garfinkle |l3| has found that the solution ( |2.10D is in good 
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agreement with the numerical results of |J at an intermediate time for the 
special value c 2 = 7/8 = 0.875, which is also close to 1. 

Prom the dependence ( |8.21| ) of the limiting horizon radius on the pertur- 
bation amplitude (proportional to the sole non-vanishing integration con- 
stant C or B), the critical exponent is 7 = 1/k. For the preferred value 
c 2 ~ 1, this leads to 7 ~ 2/5 |l4| in case a) with n = 1, or 7 ~ 1/3 in case 
b). 



9 Discussion 

The reason for writing this paper has been twofold. In the first part we have 



extended previous works |p~3| , 14] motivated by the quest for the critical so- 
lution underlying the collapse of a massless and minimally coupled scalar 
field in 2+1 dimensional anti-de Sitter spacetime. In particular, we have 
constructed a family of exact solutions with CSS behavior near the singu- 
larity and AdS behavior at spatial infinity. This one-parameter (c) family 
of time-dependent solutions interpolates between the static vacuum solution 
(c 2 = 0) and the massive BTZ black hole (c 2 = 00), with the critical value 
c 2 = 1 lying at the boundary between two different singularity types (point 
singularity, and line singularity). 

These solutions served as the starting point for the subsequent linear 
perturbation analysis. A crucial question is that of the choice of the correct 
boundary conditions to be imposed along the null boundaries u = 0, v = 0, 
and the AdS boundary x = x\. These boundary conditions have been 
"tested" in the simple case c 2 = 0, where the scalar field decouples and the 
linear perturbations describe exactly the small mass static BTZ black hole. 
In this respect, the essential difference with our previous work [|14| is that 
we now recognize that the condition (C) that the perturbations do not blow 
up at the AdS boundary is always satisfied up to gauge transformations; we 
can then impose the natural smooth matching conditions (Bl) (Eq. ( |6.8D ) 
on the boundary v = 0, instead of the weaker condition f(0) = used in |L4|] 
(resulting in significantly different behaviors of the coordinate singularity 
and apparent horizon). We also considered a new and crucial condition (D), 
stating the "regularity at birth" of the apparent horizon. Fulfilment of the 
selected boundary conditions then gives rise to two distinct growing modes 
a) (Eq. (|8.8j ) with n = 1) and b) (Eq. (^1])). The application of condition 
(D) is very suggestive because, together with the other boundary conditions 
it selects, for the mode a), a unique value for both the eigenvalue k and the 
parameter c. 
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The ensuing physical picture, along with the ambiguities inherent to 
the linear perturbations approximation, is also discussed. For the mode 
a) the singularity and the apparent horizon appear simultaneously on the 
null line it = at the time v = and evolve in the region v > in a 
physically meaningful way. On the other hand, for the mode b) (which is 
the extension to c 2 / of the single BTZ mode found for c 2 = 0), the 
singularity still appears for v = 0, but the apparent horizon seems to be 
eternal, as in the case of the static BTZ black hole. This formal solution 
cannot describe actual gravitational collapse with regular initial conditions. 
Indeed, we conjecture [22] that the solution b) is actually a linearization of 
the singular static solution (4.4). We would like to argue that ultimately the 
boundary conditions (apart from natural requirements on the consistency of 
the linear approximation) should be chosen so as to optimize the agreement 
of the perturbative picture with realistic gravitational collapse. On such 
a basis the mode b) should be excluded as unphysical. This would mean 
that there is exactly one growing mode a), showing that the background 
quasi-CSS solution is indeed the critical solution for gravitational collapse. 

An objection against this view is that the quasi-CSS solutions presented 
in this paper do not have a regular timelike origin, while all near-critical col- 
lapse scenarios do evolve about a regular timelike origin. However we argue 
that a critical solution, which lies in the boundary surface separating differ- 
ent collapsing regimes, does not necessarily share with the generic collapsing 
configurations a regular central behavior, i.e. the critical solution can lie at 
the boundary of the set of generic collapsing solutions without belonging 
to the set. Indeed, further investigations currently under progress seem to 
show that certain perturbations of our quasi-CSS solutions do present a 
regular timelike center, and that, moreover, a subset of these perturbations 
correspond precisely to extended (to A < 0) Garfinkle solutions. This would 
suggest that our quasi-CSS solutions, with singular center, could be exact 
critical solutions, with near-critical solutions well described by the extended 
Garfinkle solutions. 

We finally turn to the calculation of the critical exponent 7. Numerically, 
Pretorius and Choptuik H found 7 in the range 1.15 < 7 < 1.25, while 
Husain and Olivier || estimated 7 ~ 0.81. Choosing the critical value 
c 2 = 1 (which is also motivated by the application of boundary condition 
(D)), we obtain for the mode a) the value 7 = 0.4, which differs significantly 
from both numerical values. We cannot at present give a reason for this 
disagreement (which recalls the disagreement between the numerical Q and 
analytical Q determinations of the critical exponent for scalar field collapse 
in 3+1 dimensions). Despite this shortcoming, we find it rather satisfactory 
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to have been able to reproduce from this solution a number of characteristic 
features of critical collapse and black hole formation. 
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Figure 1: Penrose diagram of the CSS solutions eq. fl4.3| ) (A = 0) a) for 
c 2 < 1, with a point singularity (dot) at u = 0, v = +oo; b) for c 2 > 1, with 
a line singularity (double line) at u = 0. 




Figure 2: Penrose diagram of the extended quasi-CSS solutions (A < 0) a) 
for c 2 < 1; b) for c 2 > 1. 
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Figure 4: Penrose diagram of the BTZ black hole (|7.8| ). The two spacelike 
central singularities (double lines) are shielded by two null horizons (the 
diagonals). 
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Figure 6: Scalar perturbations with k = c 2 + 2 (case b)). The apparent 
horizon prexists the spacelike singularity. 
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